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ABSTRACT 

Motivated by recent observational results that focus on high redshift black holes, 
we explore the effect of scatter and observational biases on the ability to recover the 
intrinsic properties of the black hole population at high redshift. We find that scatter 
and selection biases can hide the intrinsic correlations between black holes and their 
hosts, with 'observable' subsamples of the whole population suggesting, on average, 
positive evolution even when the underlying population is characterized by no- or 
negative evolution. We create theoretical mass functions of black holes convolving the 
mass function of dark matter halos with standard relationships linking black holes 
with their hosts. Under these assumptions, we find that the local AIbh — cr correlation 
is unable to fit the z = 6 black hole mass function proposed by Willott et al. (2010), 
overestimating the number density of all but the most massive black holes. Positive 
evolution or including scatter in the Mbh — c correlation makes the discrepancy worse, 
as it further increases the number density of observable black holes. We notice that 
if the Mbh — cr correlation at z = 6 is steeper than today, then the mass function 
becomes shallower. This helps reproducing the mass function of z = 6 black holes 
proposed by Willott et al. (2010). Alternatively, it is possible that very few halos (of 
order 1/1000) host an active massive black hole at z = 6, or that most AGN are 
obscured, hindering their detection in optical surveys. Current measurements of the 
high redshift black hole mass function might be underestimating the density of low 
mass black holes if the active fraction or luminosity are a function of host or black 
hole mass. Finally, we discuss physical scenarios that can possibly lead to a steeper 
A/bh ~ cr relation at high redshift. 
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1 INTRODUCTION 



The constraints on black hole masses at the highest red- 
shifts currently probed, z ~ 6, are few, and seem to provide 
conflicting results, (i) There seems to be little or no corre- 
lation between black hole mass and velocity dispersion, a 
( [Wang et al.||20fo[ ) in the brightest radio-selected quasars, 
(ii) typically black holes are 'over-massive' at fixed galaxy 
mass/velocity dispersion compared to their z = counter- 
parts (e.g., Walter et al. 2004; at lower redshift see also 
McLure fc Dunlop|[200^ 



2006 



DecarH et al. 2010 



Shields et al.| 



Merloni 



2010 



20061 |Peng et al.| 



Woo et al. 



20081, 



but (iii) analysis of the black hole mass/luminosity function 
and clustering suggests that either many massive galaxies do 



not have black holes, or these black holes are less massive 
than expected ( [Willott et al.|2010[ WIO hereafter). 



As a result of point (ii), most authors propose that 
there is positive evolution in the Mbh— galaxy relationships, 
and quantify it as a change in normalization, in the sense 
that at fixed galaxy properties (e.g. velocity dispersion, stel- 
lar mass), black holes at high redshift are more massive 
than today. For instance, Merloni et al. (2010) prop ose that 
Mbh - M, evolves with redshift as (1 -I- z)^'^^ while |Deca"rii 



et al. (20101 suggest (l-l-z)"'^**. Point (iii) above, however, is 
inconsistent with this suggestion unless only about 1/100 of 
galaxies with stellar mass ~ 10^° — 10^^ Mq at z — Q host a 
black hole (WIO). These galaxies are nonetheless presumed 
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to be the progenitors of today's massive ellipticals, which 
typically host central massive black holes. 

When inferences on the population of massive black 
holes at the highest redshift are made, we have to take into 
consideration two important selection effects (see Lauer et 
al. 2007b). First, only the most massive black holes, power- 
ing the most luminous quasars, can be picked up at such high 
redshifts ( |Shen et al.|2008||Vestergaard et al.|2008[ ). Second, 
as a result of the limited survey area of current imaging 
campaigns, only black holes that reside in relatively com- 
mon galaxies can be recovered. Taken together, these biases 
imply that the observable population of black holes at high 
redshift will span a narrow range of masses and host prop- 
erties (see also Adelberger et al. 2005, Fine et al. 2006). 

In this paper, we explore the impact of these observa- 
tional biases on attempts to recover the intrinsic proper- 
ties of the black hole population. Our calculations are based 
on simple models grounded on empirical relations measured 
at much lower redshift, and therefore our results should be 
treated with caution. The aim of this paper is only to high- 
light the effects of the different factors that can influence 
the measurement of the intrinsic properties of the black hole 
population at high redshift. 

In section 2 we describe how we generate Monte Carlo 
realizations of the AIbh — o relation at z = 6 varying the 
slope and normalization. We then select 'observable' systems 
from these samples, considering both 'shallow' or 'pencil 
beam' surveys, and test how well we can recover the param- 
eters of the Mbh — cr relation from the 'observable' systems. 
In section 3 we discuss theoretical mass functions of black 
holes derived from the mass function of dark matter halos 
and various assumptions for the Afan — o" relationship. Using 
these results, we test what assumptions can reproduce the 
black hole mass function derived by WIO. We also discuss 
(section 4) why obtaining constraints on the average accre- 
tion rates and active fraction of black holes as a function of 
host mass is crucial to our understanding of the high-redshift 
massive black hole population. Finally, in section 5 we pro- 
pose a simple theoretical framework that leads to selective 
accretion onto black holes in a manner that reconciles the 
observational results (i)-(ii) and (iii) above. 



2 SCATTER AND EVOLUTION OF THE 
A/bh - (7 RELATION AT HIGH REDSHIFT 

We can qualitatively show the effects of selection biases with 
a simple exercise. Let us assume an evolution of the Mbh — cr 
relationship of the form: 



Mbh., 



10'' M, 



(2001^)" (^ + ^)' 



(1) 



.200kms- 

where a is a function of redshift. Let us now also assume 
that at fixed a the logarithmic scatter in black hole mass is 
A =0.25-0.5 dex {Mbh = Mbh,^ x 10^*, where S is normally 
distributed, see, e.g., Gultekin et al. 2009, Merloni et al. 
2010. The results are qualitatively unchanged for a uniform 
distribution in log A.) 

We create a Monte Carlo simulation of the A/bh — cr re- 
lation at 2 = 6 assuming different values of a and 7. For this 
exercise we run a number of realizations N{Mh) oc l/n{Mh), 
where n is the number density of halos of a given mass 
{Mh) calculated through the Press & Schechter formalism. 



We then select only systems that are likely to be observed. 
We consider a shallow survey and a pencil beam survey. 
A wide, shallow survey would preferentially select systems 
with high luminosity, but has the advantage of a large area. 
For instance the SDSS quasar catalogue selects sources with 
luminosities larger than Mi = —22.0 (~ 10*^ erg s~^) over 
an area of 9380 deg^ , corresponding to a volume of almost 
7 comoving Gpc'' at z=6. To simulate a shallow survey, we 
select black holes with a sizeable mass, implying that large 
luminosities can be achieved, Mbh > 3 x 10** Mq (see, e.g. 
Salviander et al. 2007; Lauer et al. 2007b, Vestergaard et al. 
2008, Shen et al. 2008, 2010 for a discussion of this bias), 
and hosted in halos with space density n > 1 Gpc~^. Pencil 
beam surveys can probe fainter systems, but at the cost of 
a smaller area, e.g. the 2 Ms Chandra Deep Fields cover a 
combined volume of ~ 10^ comoving Mpc'^ at 2: = 6 and 



reach flux limits of ~ 10 and 



10" 



'^^ erg cm ^ s ^ m 



the 0.5-2.0 and 2-8 keV bands, respectively (the flux limit 
corresponds to a luminosity ~ 10*"* and ~ 10^* erg s~^at 
z = 6). As an example of a pencil beam survey, we select 
black holes with mass A/bh > lO'^ Mq hosted in halos with 
density n > 10^ Gpc"^ 

To select sources that are observable in current sur- 
veys, we link the velocity dispersion, a, to the mass of the 
host dark matter halo. Empirical correlations have been 
found between the central stellar velocity dispersion and 
the asymptotic circular velocity (14) of galaxies (Ferrarese 
2002; Baes et al. 2003; Pizzella et al. 2005). Some of these 
relationships (Ferrarese, Baes) mimic closely the simple 
a — Vc/\/i deflnition that one derives assuming complete 
orbital isotropy. Indeed, it is difficult to imagine that the 
ratio between a and Vc for massive, stable systems evolves 
strongly with redshift and that it can be much different 
from \/3 or \/2 (see Binney & Tremaine 2008) . Since the 
asymptotic circular velocity (14) of galaxies is a measure of 
the total mass of the dark matter halo of the host galax- 
ies, we can derive relationships between black hole and dark 
matter halo mass, adopting, for instance. Equation 1 with 
a = A and 7 = 0: 



Afh 



:8.2xlO^^A/p 



AfB 



IQ^M, 







i irj 



^4, 187r2 



-1/2 



{1 + 2 



-3/2 



(2) 

In the above relationship, Ac is the over-density at viri- 
alization relative to the critical density. For a WMAP5 
cosmology we adopt here the fltting formula Ac = 
187r^ -f 82d - 39d^ (Bryan & Norman 1998), where d = 
Q^ — 1 is evaluated at the redshift of interest, so that 
n^ = n^(l + zf/{n,n{l + zf -h ^A -^ nfe(l + zf). Given 
the mass of a host halo, we estimate the number density from 
the Press & Schechter formalism (Sheth & Tormen 1999). 
In this section we assume a = 14 //S, where 14 is the virial 
circular velocity of the host halo. The results of this experi- 
ment are not strongly dependent on this speciflc assumption; 
in Section 3 below we discuss different scalings. Kormendy 
et al. (2011) question a correlation between black holes and 
dark matter halos (but see Volonteri et al. 2011 for an up- 
dated analysis). We notice that in any case Kormendy 's ar- 
gument is not a concern here, as at large masses Kormendy 
et al (2011b) suggest that a cosmic conspiracy causes a and 
14 to correlate, thus making the link between M and 14 ad- 
equate. In the Monte Carlo simulation, at fixed halo mass 
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Figure 1. Top panel: Mbh — f relation at z = 6, assuming a = 6, 
7 = 0, and a scatter of 0.25 dex. Cyan dots: 'observable' popula- 
tion in a shallow survey. Blue line: linear fit to this 'observable' 
population, yielding a = 2.7. Green dots: 'observable' population 
in a pencil-beam survey. Dark green line: linear fit to this 'ob- 
servable' population, yielding a = 4.5. Red line: fit to the whole 
population, yielding a = 6. Yellow dashed line: Mbh — ct at z = 
(Equation 1 with a = 4 and 7 = 0). Bottom panel: same for 
o = 8, 7 = —1. 



{Mh, hence, a), we derive Mbh, a from the adopted Mbh — cr 
relation (i.e. depending on the choice of a and 7), and then 
we draw the black hole mass from Mbh ~ Mbh.o- x 10 
with varying values of the scatter A. 

First, we test a no-evolution case, where we set a = 4 
and 7 = 0. We fit, in log-log space, the Mbh — cr relation 
of black holes implied by the 'observable' population, con- 
sidering both a shallow and pencil beam survey. In the no 
evolution case, we find Oflt — 1 and "fat — 0.7 in the 'shal- 
low' survey, with almost all 'observable' black holes lying 
above the a — 4 and 7 = line, suggesting 'overmassive' 
black holes, only because of the mass threshold that was 
imposed on the sample. In the 'pencil beam' survey we find 
Qfit — 2.5 and yat ~ —0.2. In either case, fitting only the 
'observable' population yields a much shallower slope than 
that characterizing the whole population. 

In Figure [T] we show a Monte Carlo simulation of the 
Mbh— o" relation at z = 6 one would find assuming A = 0.25, 
Q = 6 and 7 = (top panel), and a = 8 and 7 = — 1 
(bottom panel). In section 3, we will show that these par- 
ticular choices of a and 7 are motivated by our attempt to 
fit the black hole mass function of WIO. In the a = 6 and 
7 = Monte Carlo simulation, we find that the best fit has 
Qfit ^ 2.7±0.2 and 7fit ~ 0.45±0.04 for the 'shallow' survey. 
The apparent normalization of the relationship therefore in- 
creases by 0.35 dex (all the blue points lie above the yellow 
line in the top panel of figure [T| . So while the underlying 
population is characterized only by a change in slope (with 
respect to the z — relationship), what would be recovered 
from the 'observable' population is a shallower slope and a 
positive evolution of the normalization (in agreement with 



point (ii) in §1). We note, additionally, that the smaller the 
range in Mbh that is probed, the more likely it is that the 
scatter A hides any correlation, likely explaining the lack of 
correlation (point (i) in §1) found by Wang et al. (2010 rj 
If we increase the level of scatter (A) the slope of the rela- 
tionships recovered from the Monte Carlo sample becomes 
progressively shallower. 

We can repeat the same exercise for, e.g., a = 8 and 
7 = — 1, and although the underlying population has a much 
steeper slope and a negative evolution of the normalization 
of the Mbh — o" relation with redshift, the 'observable' pop- 
ulation in the shallow survey would nevertheless display no 
evolution at all (blue vs yellow lines in Fig. [l]). 

Summarizing we find that, (1) selection effects can 
severely alter the mapping between black mass and host 
galaxy velocity dispersion, leading to observed black hole 
populations that are more massive than the true distribu- 
tion. (2) Scatter and selection effects can mask correlations 
between black mass and host galaxy properties, leading to 
observed Mbh — cr relations that are shallower than the true 
relation. Although the quantitative results must be taken 
with caution, the existence of biases towards measuring a 
positive evolution in the black hole-host correlations in- 
duced by selection and scatter is generically a robust result 
(e.g.. Shields et al. 2006, Salviander et al. 2007; Lauer et al. 
2007b). 



3 IMPACT OF EVOLUTION OF Mbu cr 
RELATION AND SCATTER ON THE 
BLACK HOLE MASS FUNCTION 

We now turn to the mass function of black holes, and how 
its shape and normalization are affected by the evolution 
of Mbh — cr relation and its scatter. We create theoretical 
mass functions based on Equation 1 coupled with the Press 
& Schechter formalism, exploring how different values of a 
and 7 influence its functional form. As discussed in section 
2, one can derive relationships between black hole and dark 
matter halo mass given a relationship between black hole 
mass and velocity dispersion (Equation 1), a relationship 
between velocity dispersion (cr) and asymptotic circular ve- 
locity (virial velocity, Vc), and the virial theorem. For in- 
stance, assuming Equation 1 with q = 4 and 7 = 0, and 
o" ~ Vc/\/3 one derives Eq. |2l while if we assume the rela- 
tionship proposed by Pizzella et al. (2005) between a and 



Vc-. 



Mh 



:4.1XlO^^Mp 



Mb 



lO^M, 







i ijl 



n^, 187r2 



-1/2 



(1+^) 



-3/2 



(3) 

To consider the range of possible black hole mass func- 
tions, we adopt the two mappings between black hole mass 
and halo mass provided by Equations [2] and |3] We first con- 
sider the resulting black hole mass function when we adopt 
the local Mbh — cr relation (a = 4 and 7 = 0), and we will 



^ Wang et al. did not attempt any fit to the Mbh — "" relation. 
They note that they find significant scatter, extending to over 
3 orders of magnitude, and that most of the quasar black hole 
masses lie above the local relationship. See also Shields et al. 
(2006) for quasars at 2 = 3. 
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then investigate how the mass function changes if we vary 
a and 7. In particular, we will focus on a = 6 and 7 = 0, 
and a — S and 7 — —1, because, as shown below, a steeper 
Mbu — o" relation yields better agreement between theoret- 
ical black hole mass functions and WIO. 

We can estimate the mass function of black holes by 
convolving equations pjis] (and their possible redshift evolu- 
tion) with the mass density of dark matter halos with mass 
Mh derived from the Press & Schechter formalism (Sheth & 
Tormen 1999): 



dN 



dN dlogMh 
d log Mbh ^ d log Mh d log Mbh ' 



(4) 



We assume for the time being that black holes exist 
in all galaxies. The effect of dropping this assumption is 
discussed in detail in Section 4. In figure p] we compare the 
mass function derived using this technique to the mass func- 
tion proposed by WIO, based on the luminosity function of 
quasars selected by the Canada-France High-z Quasar Sur- 
vey, assuming a duty cycle (corresponding to the fraction of 
black holes that are active, we will refer to this quantity as 
the active fraction, AF, below) of 0.75 and assuming a log- 
normal distribution of Eddington fractions, /Edd, centered 



at 0.6 with standard deviation of 0.30 dex (see also Shankar 
et al.|[2010a K WIG further assume the same fraction of ob- 
scured AGN as observed at lower redshift (z = — 2, Ueda 
et al. 2003), and correct for Compton thick AGN following 
Shankar et al. (2009) . Note that the evolution of the frac- 
tion of obscured or Compton thick AGN (currently not well 
constrained at high redshift, but see Treister et al. 2011) can 
strongly influence the results by hiding part of the black hole 
population (see section 4). 

If the Mbh — o" relation evolved with redshift as pro- 
posed by Woo et al. (20081, 7 = 3.1, the number density 
of black holes in the mass range 10^ — 10^ M0 would be 
~ 0.5 and 10~* comoving Mpc~^ respectively (the curve 
corresponding to this very strong evolution is not shown in 
the figure). We note, however, that the sample analyzed by 
Woo et al. (2008) is at z ~ 0.4, and there is no guarantee 
that such evolution holds at higher redshift. 

In all cases the analytical models greatly over-estimate 
the mass function at masses Mbh < 10® Mq, and possibly 
at all masses - when we add the suggested positive redshift 
evolution of the Mbh— galaxy relationships. 

In figure [3] we show instead the mass function we find 
when we assume different a and 7 values, with and without 
scattering. We include scattering, at the level of A = 0.5, by 
performing a Monte Carlo simulation, where for each black 
hole mass we create 500 realizations of the host mass. The 
WIO black hole mass function can be reproduced by a simple 
model that has a = 8 and 7 = —1, if no or little scatter in 
the black hole properties with galaxy mass is present. We 
see that as a increases the mass function becomes shallower. 
At fixed black hole mass, above the 'hinge' of Equation 1 
(200 kms"'^) black holes will be found in comparatively less 
massive galaxies, that have a higher density. On the other 
hand, below the 'hinge', the host of a black hole of a given 
mass would be more massive than in in the a = 4 case, 
hence with a lower space density. This effect makes the mass 
function shallower. Any decrease in 7 tends to shift the black 
hole mass function to lower number densities at all masses. 

However, a significant amount of scatter increases the 
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Figure 2. Mass function of black holes. Black dots: Willott et al 
2010. Orange stars: Mbh - M, + Stark et al. 2009 (see Willott et 
al. 2010 for details). Blue long dashed curve: Press & Schechter 
-I- equation [2] (a = 4). Dark green short dashed curve: Press & 
Schechter -|- equation [3](o = 4). 



number density of observable black holes, as shown in the 
bottom panel of Figure p] This effect has been discussed 



extensively by Lauer et al. (20071, and we refer the reader 



to this paper for an exhaustive demonstration of its conse- 
quences. iLaueF^F^L] ((2007b start from the luminosity func- 
tion of galaxies, rather than the meiss function of dark matter 
halos, and the fact that the most luminous galaxies are in the 
exponential part of the luminosity function, implies that the 
scattering of very high-mass black holes (Mbh — 10® Mq) in 
lower mass galaxies has a stronger effect than the scattering 
of low mass black holes in larger galaxies. A similar con- 
clusion applies to the mass function of dark matter halos. 
Additionally, since the halo mass function becomes expo- 
nential at lower masses at high redshift, the effect of scatter 
on the shape of the black hole mass function becomes no- 
ticeable already at A/bh — 10^ Mq. Including scatter, the 
simple model with a — 8 and 7 = — 1 is now a much poorer 
fit to WIO mass function, but it still reproduces their slope 
very well. 

Summarizing, we find that the local Mbh — o" relation 
{a — 4 and 7 = 0) is unable to reproduce WIO results, 
even more so when a level of scatter compatible with ob- 
servational results (A — 0.25 — 0.5) is included. A steeper 
-Ai^BH — o" relation, with possibly a negative evolution (e.g., 
Q = 8 and 7 = — 1) provides a better fit, although high 
levels of scatter require an even more dramatic steepening 
of the slope in order to match the mass function proposed 
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Figure 3. Mass function of black holes. Top panel: we vary the 
slope (a) and normalization (7) of the Mbh — cr relation, and 
assume no scatter in the relationship (A = 0). Bottom panel: we 
vary the slope of the Mbh — o" relation, and include scatter in the 
relationship (A = 0.25 dex or A = 0.5 dex). Black dots: Willott 
et al 2010. All curves assume Press & Schechter + equation [2] 
with varying a as labelled in the Figure. 



by WIO. While the direct comparison with WIO strongly 
depends on the limitations of our empirical model, the re- 
lationship between increased scatter in the Mbh — o" and 
increased number density of black holes is a robust result 
that directly follows from the analysis presented in |Lauer| 



4 OCCUPATION FRACTION OF QUIESCENT 
AND ACTIVE BLACK HOLES 

In section 3 we demonstrated how we derive the theoreti- 
cal mass function of black holes from the mass function of 
their host halos and the relation between black hole and 



halo masses (e.g., Haiman & Loeb 1998 Wyithe & Loeb 
2002). However, when we convolve equations [2] and |3| with 
the mass density of dark matter halos to derive a black hole 
mass function, we have to make a conjecture about the frac- 
tion of halos of a given mass which host a black hole, the 
occupation fraction (OF): 

d log AfBH " ^ "' ^' d log A4 d log Mbh ' ^^ 

In the top panel of Fig.|4]we show black hole mass func- 
tions resulting from different choices of the OF. It is clear 



that a decreasing OF can compensate for an increased scat- 
ter in shaping the black hole mass function. As a result of 
this degeneracy, we can reproduce the WIO mass function 
for a range of values for the OF and scatter. Even adopt- 
ing a = 4 (slope of Mbh — cr relation at the present day) 
with a sensible scatter can fit the data, at the cost, how- 
ever, of making the presence of a black hole (regardless of 
its shining as a quasar) a very rare instance. We note that 
it is conceivable that the OF is not constant over all host 
masses, and a non trivially constant OF is expected par- 
ticularly at high-redshift, close to the epoch of galaxy and 
black hole formation ( Menou et al.|2001 1. At face value, the 
WIO data can be reproduced by OF= Mh/5 x 10^^ Mofor 
Q = 4, 7 = and A = 0.25, or OF= (Mh/IO^MQ)^'^^ for 
Q = 4, 7 = and A = 0. Such occupation fractions are 
several orders of magnitude lowers than predicted by mod- 
els of formation and cosmic evolution of black holes (e.g., 
Volonteri||2010l, and we therefore still prefer solutions with 



steeper slopes. We will explore self-consistently OF and its 
relationship with the establishment of the A^bh — o" relation 
as a function of black hole formation and growth physics in 
a future paper. 

Throughout this paper we have compared our theoreti- 
cal mass function of black holes to constraints derived indi- 
rectly from the luminosity function of quasars in WIO, rather 
than from direct black hole mass measurements (that are 
rather unfeasible at 2: = 6). Empirically, one can derive the 
mass function of black holes from the luminosity function of 
quasars and a relation between black hole mass and quasar 
luminosity (e.g., |Shankar et al.|2010a|b| [Willott et al.|20"l0l ): 



dN 



dN dlogL 



d log A/bh d log L d log A/bh 



(6) 



For instance, we can estimate the mass function of black 
holes from the bolomotric luminosity of radio-quiet quasars 
(Hopkins et al. 2007) assuming (1) that all black holes are 



active, (2) that all black holes radiate at the same Edding- 
ton fraction, /Bdd (based on various observational results we 
expect high redshift quasars to radiate close to the Edding- 
ton limit, see WIO and references therein). The mass of a 
black hole powering a quasar with luminosity L is then: 



Mbh 
lO^Mp 



3 X 10" 



./Edd ig 



(7) 



and one can trivially turn the luminosity function into 
a mass function. As discussed by WIO, their mass function 
is derived assuming similar values of the Eddington ratio 
and the active fraction using a more accurate technique (see 
[Shankar et al.|2010a[ for details) . Our simple approach pro- 
vides results consistent with WIO if we assume a constant 
/Edd ~ 1. 

When we deconvolve the luminosity function of quasars 
to derive the black hole mass function we have to assume an 
active fraction (AF): 



(8) 



dN dN dlogL 

d log Mbh " ^ ' ""' dAogL d log Afen ' 

where we indicate that both the active fraction and the 
Eddington ratio can be functions of the black hole (and host) 
properties, and of cosmic time. 
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The intrinsic shape of the mass function changes as a 
function of AF and /Edd, and, in particular, any departure 
from the assumptions /Edd ~ 1 and AF= 1 (that are up- 
per limits to both quantities) will drive the mass function 
of black holes 'up', that is, they will increase the number 
of black holes at a given mass. We have therefore to bear 
in mind that the semi-empirical mass function derived by 
WIO might be underestimating the mass function. The lower 
panel of Fig. |4] shows how the mass functions one derives 
from a luminosity function are modified by an AF or /Edd 
that depend on the BH mass. For instance, we can trivially 
assume that /Edd = Mbh/W^ Mq for Mbh < 10* M©, and 
/Edd = 1 otherwise. Then Mbh = lO*M0(L/3 x 10^^1/0)° ■^ 
using equationlTlin the last expression. In the same figure we 
also show the effect of a mass dependent AF, where we adopt 
the simple expression AF= Mbh/ 10^ Mq for Mbh < 10* 
Mq, and AF= 1 otherwise. These specific forms of the mass 
dependence of AF and /Edd are motivated by the expecta- 
tion that the most massive black holes at the earliest cos- 
mic times are all actively and almost constantly accreting 



| Haiman|2004l|Shapirop005l|Volonteri fc Rees|2005l [2006 1 
Such functional forms are here only used to prove that non 
constant accretion rates modify the expectations in terms 
of mass function and active fraction, but the expressions we 
adopt should be considered representative of any class of 
accretion rates and active fractions that are not constant, 
rather than actual predictions. 

If the Eddington ratio and/or AF are a function of halo 
or black hole mass, then what one derives from flux-limited 
surveys will be dependent on a combination of various prop- 
erties. FigurelSlshows simple examples. We build a sample of 
black holes and hosts by performing a Monte Carlo sampling 
as described in Section 3 for two A/bh — cr relations (a = 4 
and a = 6, both with 7 = 0) each with a scatter of 0.25 
dex. We assign to each black hole a luminosity by assuming 
either a constant /Edd = 1, or that /Edd ~ Mbh/W^ Mq or 
Mbh < 10^ Mq and /Edd =1 for Mbh > 10* M0(we note 
that the assumption that /Edd scales with the halo mass, 
rather than the black hole mass, yields very similar results). 
Even for a constant Eddington ratio, scatter in the A/bh — o" 
relation implies that at fixed halo mass the black hole mass 
and hence its luminosity is not univocally determined. The 
'observed' active population of black holes will therefore be 
different from the 'intrinsic' active fraction, which in this ex- 
ercise was set to unity. A comparison between the left and 
right panels underscores how the Mbh — cr relation itself 
shapes the fraction of active black holes which are detected 
in optical imaging surveys. 

We notice that obscuration plays a role similar to the 
occupation or active fraction. If, say, all black holes are ac- 
tive, but a large fraction are Compton thick, then a large 
population of obscured quasars would be unaccounted for in 
optical quasar surveys. There is indeed evidence for a large 
fraction of high rcdshift quasars being obscured (e.g. 
[Franca et aL[2005^ Treister et al.,2009|[20TT| ). 



5 ACCRETION EFFICIENCY AND HOST 

MASS 

In the previous section we discussed how the 'observed' frac- 
tion of active black holes, which goes into determining the 
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Figure 4. Top: theoretical mass function of black holes derived 
from the mass function of dark matter halos (Eq. 5). Black dots: 
Willott et al 2010. Other curves as marked in the figure (from top 
to bottom). Willott et al. results can be reproduces for a range 
of possible assumptions on the relation between holes and halos. 
Bottom: empirical mass function of black holes derived from the 
luminosity function of quasars (Eq. 7 and Eq. 8, adopting the 
bolometric luminosity function of Hopkins et al. 2007). Triangles: 
/Edd = 1 a^nd AF = 1 . The blue solid and red dashed curves show 
how mass-dependent luminosities or active fractions can modify 
the shape of the mass function one would infer. 



'observed' mass function depends on the Mbh — o" relation- 
ship and on the link between accretion rate and black hole- 
host masses. In this section, we explore the consequences 
and likelihood of a galaxy mass-dependent black hole ac- 
cretion rate. This hypothesis is plausible, as the gas supply, 
especially at high redshift is likely dependent on the envi- 
ronment and mass of the host. For instance cold gas that 
flows rapidly to the center of galaxies from filaments around 
haloes plays a major role in the buildup of massive galaxies 
at high redshift ( [Brooks et al.|2009|[Governato et al.|2009[ ), 
with a transition expected to occur when a galaxy has mass 
above lO^^Af0, where gas is shocked before it can reach the 
galaxy's disk. Cold gas fiowing into halos along large-scale 
structure filaments may however be dense enough to pen- 
etrate the shock front and deliver cold gas to the galaxy. 
Galaxies that form within a gas-rich filament will accrete 
gas from this cold flow and grow substantially before the fil- 
ament dissipates. These galaxies embedded in filaments are 
expected to be high peaks of the density field, hence among 
the most massive at early times. 

Additionally, as discussed in section 3, instead of an 
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Figure 5. Fraction of active blacl< holes than can be detected in 
a survey with given luminosity and volume limit. Red histograms 
(leftmost side of each panel): black holes with Lmin = lO** erg/s 
and riniin = 10~^ Mpc""^ (example of a pencil beam survey). 
Black histograms (rightmost side of each panel) : black holes with 
^min = lO'l'' erg/s and riniin = 10~® Mpc~^ (example of a shallow 
survey). Left panels: o = 4; A =0.25 dex. Right panels: a = 6; 
A =0.25 dex. Top panels: /Edd = 1- Bottom panels: /sdd = 
Mbh/10» Mq. 



overall normalization evolution, the link between black holes 
and their hosts might be better explained by an evolution in 
the slope of the Mbh — o" relationship. We are not claiming 
here that the evolution has the exact form that we use in 
this paper (Equation [l|. We here discuss a possible physical 
scenario that can lead to a steeper Mbh — cr relation. 

In the following toy model we just explore what physical 
process could drive the establishment of a given Mbh — cr 
relation at a given redshift {z — 6 in this particular case). 
In other words, if the slope of the Mbh — o" relation at 2 = 6 
has a given a, what can be the driver of such correlation? 

Let us assume that all black holes start with the same 
initial mass. Mo, and let them grow until z — 6 {Ih = 0.9 
Gyr). If at z = 6 the slope of the Mbh — cr relationship is 
Q, and we assume that on average black holes accrete at a 
fraction (/Edd) of the Eddington rate, then we can relate this 
average accretion rate of the mass of the hole, Mbh, 6 and 
the velocity dispersion, ae, of the host at 2 = 6, as follows: 



Mbh,6 = Mo exp ( (/Edd) 
where tEdd — Afenc^/I/Edd = 



tH 1 

^Edd 



(9) 



. r-— = 0.45 Gyr (c is 

4-77 G rrip •> ^ 

the speed of light, or is the Thomson cross section, m^ is 



the proton mass), and the radiative efficiency, e ~ 0.1. We 
can also express the relationship between black hole mass 
Mbh, 6 and ae at 2 = 6 as: 



0-6 = 200 km s 



-1 ( Mb 



VIO8M0 



1/0 



(10) 



so that the average accretion rate for a black hole that 
grows within a galaxy that has a velocity dispersion cts at 
2 = 6 results; 



(/Edd) 



^Edd 



tH 1 



In 



Afo 



0-6 



200 km s- 



(11) 



Equations 9-11 are based on the initial conditions and on 
the properties at 2 = 6 only. The accretion rate is in princi- 
ple galaxy-mass dependent, and specifically dependent also 
on the mass growth of the host. However, for the sake of 
simplicity it is here set to the average over the integration 
time. 

The average accretion rate of Eq. 11 is shown in Fig- 
ure |6] for a = 4, a = 6 and = 8. Figure [6] implies that 
if only black holes in hosts above a certain velocity disper- 
sion, or mass, or depth of the potential well, can accrete ef- 
ficiently, it is only natural to expect a different slope of the 
Mbh — <y relationship in dependence of the exact threshold. 
One possibility is that black hole growth is indeed ineffi- 
cient in low-mass galaxies at early cosmic times, because of 
the fragile environment where feedback can be very destruc- 



tive ([Milosavljevic et "ar||2009[ [Alvarez et aL]|2009[ [Park fc 



Ricotti 2010 Johnson et al. 20101. We note that all these 



possible scalings of accretion rate with halo mass are con- 
sistent with the independence of /Edd on luminosity (and 
black hole mass) found by WIO, as all their black holes have 
masses above 10* M0 , where it is indeed expected that the 
accretion rate can reach /Edd — 1 as one can infer from 
Figure [6] 

This exercise is not meant to suggest that the typical 
accretion rate has the exact value of Eq. 9, but that «/ accre- 
tion is more efficient in more massive halos, then a increases, 
while if the accretion rate is mass independent, e.g. is con- 
stant in all hosts, then a tends to lower values. 

Equation [TT] simply demonstrates mathematically that 
in order to achieve a very steep Mbh — o" accretion in massive 
haloes has to be more efficient than in small haloes ('selec- 
tive accretion'), and it should not be used to make predic- 
tions about the accretion/growth history of black holes. To 
test this suggestion, dedicated simulations that can resolve 
the growth of black holes in cosmological simulations as a 
function of the host mass are required. However, simulations 
that explore the cosmic evolution of accretion efficiency, tak- 
ing into consideration feeding and feedback, as a function of 
host mass at sufficient resolution are not currently available. 

The experiment t hat is closest in spirit to what we pro- 
pose was performed by Pelupessy et al.[ ([2007[p]who suggest 
that the more massive the host halo, the higher the Edding- 
ton fraction. A very simple fit from their simulation results 
at 2 = 6 (Figure 7 in Pelupessy et al. 2007) gives: 



■^ Indirectly, similar information can be extracted from |Sijacki| 
jet al.| ( [2009| l, although all their information on the accretion rate 
is cast in terms if black hole masses, rather than host properties. 
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Mh 



10" M 







(12) 



This equation is shown in Figure [6] for different hosts, 
where we assumed a = Vc//3, and that the Eddington limit 
is capped at /Edd ~ 1 (in the bottom panel the black hole 
mass uses Equation S|. This scaling leads to very steep re- 
lationships between black holes and hosts, as shown in the 
bottom panel of Figure[6l as the black hole mass depends ex- 
ponentially on the halo mass (if we insert equation 12 into 
equation 9) or on the time evolution of the halo mass (if 
we insert equation 12 into the expression for the accretion 
rate in Eddington units: M = /Edd(Af/iEdd)[(l — e)/e] ~ 
(Mh(t)/10")(M/tBdd)[(l - e)/e]. To integrate this equation 
properly one needs to know the growth rate of the host dark 
matter halo mass as a function of time in a ACDM cosmol- 
ogy, Mh(t). Such exercise requires either merger trees that 
track the cosmic history of dark matter halos, or an analyt- 
ical fit to their growth histories, and it is beyond the scope 
of this paper). Selective accretion, modulated by the host's 
potential and environment, is a possible key to explaining a 
shallow high redshift black hole mass function without re- 
quiring an unrealistically low occupation fraction. 

At late cosmic times we expect the interaction of black 
holes and galaxies to become more closely linked to bary- 
onic processes (e.g., bulge formation) rather than being re- 
lated to the halo mass. For instance, secular effects might 
at late times decouple the properties of the central stellar- 
dominated region from the overall dark matter halo. As an- 
other example, gas accretion through cold fiows in filaments 
is expected to occur only at early times. We can, for instance, 
see a parallel between the black hole-halo relationship and 
the baryon-halo relationship. It is expected that at very high 
redshift most halos possess a baryon fraction of order the 
cosmic baryon fraction, while at later times the baryonic 
content evolves under the effect of baryonic physics. In the 
same way at late times we expect black hole growth to be 
more closely related to the baryonic content of a galaxy, and 
hence less "selectively" linked to the host halo (cf. Volonteri 
et al. 2011). 



6 SUMMARY AND CONCLUSIONS 

Recent observational results that focus on high-redshift 
black holes provide seemingly conflicting results. In partic- 
ular: 

(i) there seems to be little or no correlation with velocity 
dispersion, a in the brightest radio-selected quasars, 

(ii) typically black holes are 'over-massive' at fixed galaxy 
mass/velocity dispersion compared to their z = counter- 
parts, 

(iii) clustering and analysis of the mass/luminosity function 
suggest that either many massive galaxies do not have black 
holes, or these black holes are less massive than expected. 

To try and understand these observational results, we 
explore the role of scatter and observational biases in recov- 
ering the intrinsic properties of the black hole population. 
We generate Monte Carlo realizations of the Mbh — o" rela- 
tion at 2 = 6, varying the slope and normalization, and select 
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Figure 6. Top: Eddington fraction as a function of halo velocity 
dispersion, derived from equation | 1 1 [ (black solid: o = 4; blue long 
dashed: a = 6; red short dashed: a = 8) and equation |12| (green 
dot-dashes). Bottom: black hole mass versus a a.t z = 6 assuming 
/Edd is a function of halo properties, and letting the holes accrete 
for tH = 0.9 Gyr. 



'observable' systems from these samples, considering either 
'shallow' or 'pencil beam' surveys. We test how well we can 
recover the parameters of the Mbh — o" relation from the 
'observable' systems only. We then create theoretical mass 
functions of black holes from the mass function of halos and 
Mbh — o" and test what assumptions can reproduce the mass 
function derived by WIO. Our techniques are very simplified 
and we use empirical correlations that are not guaranteed to 
hold at all masses and reshifts. Therefore one should not in- 
terpret our results as solutions to the three conflicting points 
mentioned above, but rather regard them as a way for under- 
standing how different physical parameters may affect black 
hole related quantities and their measurements. 
Our results can be summarized as follows: 

• Scatter and bias selections can hide the intrinsic cor- 
relations between holes and hosts. When selecting within a 
small range in black hole and galaxy masses, at the high- 
mass end, the scatter washes out correlations (see point (i) 
above), and most of the selected systems will tend to lie 
above the underlying correlation. The correlations recovered 
from 'observable' sub-samples of the whole population can 
therefore suggest positive evolution even when the underly- 
ing population is characterized by no or negative evolution. 

• The slope and normalization of the local Mbh — o" corre- 
lation are unable to produce a black hole mass function com- 
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patible with WIO, as the theoretical mass function greatly 
overstimates the density of black holes with Mbh < 10^. 
The discrepancy can be minimized if very few halos (of or- 
der 1/1000) host a massive black hole or an AGN at 2 = 6, 
or most AGN at these redshift are obscured. Mq. 

• If the Mbh — o" correlation were steeper at z = 6, then at 
fixed black hole mass high mass black holes would reside in 
comparatively less massive galaxies than in the a = 4 case. 
Their number density is therefore increased. Viceversa, low 
mass black holes would be hosted in comparatively larger 
galaxies (compare red and yellow lines in Figure 1) with a 
lower space density. This effect helps reproducing the mass 
function of 2 = 6 black holes proposed by WIG. 

• On the other hand, scatter in the Mbh — cr, at the level 
of what is observed locally, exacerbates the discrepancy, as 
it increases the number density of black holes at A/bh > 
lO'^ Mq. Any type of positive evolution of the Mbh — o" 
exacerbates this discrepancy. 

• Analysis of AGN samples might be underestimating the 
black holes mass function at low masses if the active fraction 
or luminosity are a function of host or black hole mass. 

In the near future the synergy of JWST and ALMA can 
zoom in on quasars and their hosts respectively informing 
us of their relationship and how the Mbh — o" relation is 
established, or how the accretion properties depend on the 
black hole or halo mass. In the near-IR, JWST will have the 
technical capabilities to detect quasars at z > 6 down to a 
mass limit as low as 10^ — 10^ Mq, owing to its large field 
of view and high sensitivity. At the expected sensitivity of 
JWST, ~ 1 n jy, almost 7x 10^ deg~ ^ sources at z > 6 should 
be detected ( Salvaterra et al!]|2007 l. At the same time, the 
exquisite angular resolution and sensitivity of ALMA can be 
used in order to explore black hole growth up to high redshift 
even in galaxies with high obscuration and active star for- 
mation. To date the best studies of the hosts of z ~ 6 quasar 
have been performed at cm- wavelength (Walter et al.|2004 



Wang et alT]|2010[ ). The best studied case is J1148-I-5251 at 
6.42. The host has been detected in thermal dust, non- 



thermal radio continuum, and CO line emission ( Bertoldi et 
aL][2003l [Carilh et al.|20"04} [Walter et al.|2004| ). ALMA will 
be able to detect the thermal emission from a source like 



J1148-I-5251 in a few seconds at sub-kpc resolution (Car- 
illi et al.' 2008). On a similar time-frame Dark-Energy ori- 



ented survey will provide an enormous amount of quasar 
data as ancillary science (e.g. DES, LSST). Coupling the in- 
formation we derive from these extremely large yet shallow 
surveys, with that derived from deep pencil beam surveys 
will undoubtedly deepen our understanding of the growth of 
high-redshift black holes. 

In the meantime, we need to develop dedicated cosmo- 
logical simulations of black hole formation and early growth 
that can aid the interpretation of these data. The sugges- 
tion that the accretion rate of massive black holes depends 
on their environment, (through the host halo and its cos- 
mic bias) must be tested with cosmological simulations that 
implement physically-motivated accretion and feedback pre- 
scriptions. We also need to derive predictions for the occupa- 
tion fraction of black holes in galaxies based on black hole 
formation models (Bellovary et al. in prep.). This will be 
a significant improvement over current simulations of black 
hole cosmic evolution that typically place black holes in ha- 



los growing above some threshold mass, typically ~ 10^ 
Mq , leading to a trivial occupation fraction function. There 
is no strong physical reason to believe that all and only 
galaxies with mass > W^^Mq host massive black holes in 
their centers. 

In this paper we focused on the very high-redshift Uni- 
verse, z ~ 6. Although this redshift range is not a special 
place, the concurrence of theoretical arguments and observa- 
tional constraints allow us to make simplifying assumptions 
that are not expected to be valid at later times. For instance, 
a timescale argument requires black holes to grow fast to 
reach the masses probed by current luminosity functions. In 
turn, this argument, coupled to current observational con- 
straints suggest that the most luminous quasars accrete close 
to Eddington, and that both active fraction and occupation 
fractions must be of order unity at the high-luminosity, high- 
mass end. This is not true at z — 2, where more variables 
enter into play, and make the analysis less constraining. An 
example of a detailed study that connects the mass function 
of black holes derived from the Press & Schechter formalism 
to that derived from the luminosity function is presented in 
Croton (2009). 
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